Metagenomic insights into the taxonomy, function, and dysbiosis of prokaryotic communities in octocorals.

BACKGROUND: In octocorals (Cnidaria Octocorallia), the functional relationship between host health and its symbiotic consortium has yet to be determined. Here, we employed comparative metagenomics to uncover the distinct functional and phylogenetic features of the microbiomes of healthy Eunicella gazella, Eunicella verrucosa, and Leptogorgia sarmentosa tissues, in contrast with the microbiomes found in seawater and sediments. We further explored how the octocoral microbiome shifts to a pathobiome state in E. gazella. RESULTS: Multivariate analyses based on 16S rRNA genes, Clusters of Orthologous Groups of proteins (COGs), Protein families (Pfams), and secondary metabolite-biosynthetic gene clusters annotated from 20 Illumina-sequenced metagenomes each revealed separate clustering of the prokaryotic communities of healthy tissue samples of the three octocoral species from those of necrotic E. gazella tissue and surrounding environments. While the healthy octocoral microbiome was distinguished by so-far uncultivated Endozoicomonadaceae, Oceanospirillales, and Alteromonadales phylotypes in all host species, a pronounced increase of Flavobacteriaceae and Alphaproteobacteria, originating from seawater, was observed in necrotic E. gazella tissue. Increased abundances of eukaryotic-like proteins, exonucleases, restriction endonucleases, CRISPR/Cas proteins, and genes encoding for heat-shock proteins, inorganic ion transport, and iron storage distinguished the prokaryotic communities of healthy octocoral tissue regardless of the host species. An increase of arginase and nitric oxide reductase genes, observed in necrotic E. gazella tissues, suggests the existence of a mechanism for suppression of nitrite oxide production by which octocoral pathogens may overcome the host's immune system. CONCLUSIONS: This is the first study to employ primer-less, shotgun metagenome sequencing to unveil the taxonomic, functional, and secondary metabolism features of prokaryotic communities in octocorals. Our analyses reveal that the octocoral microbiome is distinct from those of the environmental surroundings, is host genus (but not species) specific, and undergoes large, complex structural changes in the transition to the dysbiotic state. Host-symbiont recognition, abiotic-stress response, micronutrient acquisition, and an antiviral defense arsenal comprising multiple restriction endonucleases, CRISPR/Cas systems, and phage lysogenization regulators are signatures of prokaryotic communities in octocorals. We argue that these features collectively contribute to the stabilization of symbiosis in the octocoral holobiont and constitute beneficial traits that can guide future studies on coral reef conservation and microbiome therapy. Video Abstract.

Description of Alicyclobacillus montanus sp. nov., a mixotrophic bacterium isolated from acidic hot springs.

Three morphologically similar thermo-acidophilic strains, USBA-GBX-501, USBA-GBX-502 and USBA-GBX-503T, were isolated from acidic thermal springs at the National Natural Park Los Nevados (Colombia). All isolates were spore-forming, Gram-stain-positive and motile, growing aerobically at 25-55 °C (optimum ~45 °C) and at pH 1.5-4.5 (optimum pH ~3.0). Phylogenetic analysis of the 16S rRNA gene sequences of these isolates showed an almost identical sequence (99.0 % similarity) and they formed a robust cluster with the closest relative Alicyclobacillus tolerans DSM 16297T with a sequence similarity of 99.0 %. Average similarity to other species of the genus Alicyclobacillus was 93.0 % and average similarity to species of the genus Effusibacillus was 90 %. In addition, the level of DNA-DNA hybridization between strain USBA-GBX-503T and Alicyclobacillus tolerans DSM 16297T was 31.7 %. The genomic DNA G+C content of strain USBA-GBX-503T was 44.6 mol%. The only menaquinone was MK-7 (100.0 %). No ω-alicyclic fatty acids were detected in strain USBA-GBX-503T, and the major cellular fatty acids were C18 : 1ω7c, anteiso-C17 : 0 and iso-C17 : 0. Based on phenotypic and chemotaxonomic characteristics, phylogenetic analysis and DNA-DNA relatedness values, along with low levels of identity at the whole genome level (ANIb and ANIm values of <67.0 and <91.0 %, respectively), it can be concluded that strain USBA-GBX-503T represents a novel species of the genus Alicyclobacillus, for which the name Alicyclobacillus montanus sp. nov. is proposed. The type strain is USBA-GBX-503T (=CMPUJ UGB U503T=CBMAI1927T).

Complete Genome Sequence of Alkaliphilus metalliredigens Strain QYMF, an Alkaliphilic and Metal-Reducing Bacterium Isolated from Borax-Contaminated Leachate Ponds.

Alkaliphilus metalliredigens strain QYMF is an anaerobic, alkaliphilic, and metal-reducing bacterium associated with phylum Firmicutes QYMF was isolated from alkaline borax leachate ponds. The genome sequence will help elucidate the role of metal-reducing microorganisms under alkaline environments, a capability that is not commonly observed in metal respiring-microorganisms.

Complete Genome Sequence of Anaeromyxobacter sp. Fw109-5, an Anaerobic, Metal-Reducing Bacterium Isolated from a Contaminated Subsurface Environment.

We report the genome sequence of Anaeromyxobacter sp. Fw109-5, isolated from nitrate- and uranium-contaminated subsurface sediment of the Oak Ridge Integrated Field-Scale Subsurface Research Challenge (IFC) site, Oak Ridge Reservation, TN. The bacterium's genome sequence will elucidate its physiological potential in subsurface sediments undergoing in situ uranium bioremediation and natural attenuation.

The genome of the obligately intracellular bacterium Ehrlichia canis reveals themes of complex membrane structure and immune evasion strategies.

Ehrlichia canis, a small obligately intracellular, tick-transmitted, gram-negative, alpha-proteobacterium, is the primary etiologic agent of globally distributed canine monocytic ehrlichiosis. Complete genome sequencing revealed that the E. canis genome consists of a single circular chromosome of 1,315,030 bp predicted to encode 925 proteins, 40 stable RNA species, 17 putative pseudogenes, and a substantial proportion of noncoding sequence (27%). Interesting genome features include a large set of proteins with transmembrane helices and/or signal sequences and a unique serine-threonine bias associated with the potential for O glycosylation that was prominent in proteins associated with pathogen-host interactions. Furthermore, two paralogous protein families associated with immune evasion were identified, one of which contains poly(G-C) tracts, suggesting that they may play a role in phase variation and facilitation of persistent infections. Genes associated with pathogen-host interactions were identified, including a small group encoding proteins (n = 12) with tandem repeats and another group encoding proteins with eukaryote-like ankyrin domains (n = 7).

Experimental determination and system level analysis of essential genes in Escherichia coli MG1655.

Defining the gene products that play an essential role in an organism's functional repertoire is vital to understanding the system level organization of living cells. We used a genetic footprinting technique for a genome-wide assessment of genes required for robust aerobic growth of Escherichia coli in rich media. We identified 620 genes as essential and 3,126 genes as dispensable for growth under these conditions. Functional context analysis of these data allows individual functional assignments to be refined. Evolutionary context analysis demonstrates a significant tendency of essential E. coli genes to be preserved throughout the bacterial kingdom. Projection of these data over metabolic subsystems reveals topologic modules with essential and evolutionarily preserved enzymes with reduced capacity for error tolerance.

Cross-talk between catalytic and regulatory elements in a DEAD motor domain is essential for SecA function.

SecA, the motor subunit of bacterial polypeptide translocase, is an RNA helicase. SecA comprises a dimerization C-terminal domain fused to an ATPase N-terminal domain containing conserved DEAD helicase motifs. We show that the N-terminal domain is organized like the motor core of DEAD proteins, encompassing two subdomains, NBD1 and IRA2. NBD1, a rigid nucleotide-binding domain, contains the minimal ATPase catalytic machinery. IRA2 binds to NBD1 and acts as an intramolecular regulator of ATP hydrolysis by controlling ADP release and optimal ATP catalysis at NBD1. IRA2 is flexible and can undergo changes in its alpha-helical content. The C-terminal domain associates with NBD1 and IRA2 and restricts IRA2 activator function. Thus, cytoplasmic SecA is maintained in the thermally stabilized ADP-bound state and unnecessary ATP hydrolysis cycles are prevented. Two DEAD family motifs in IRA2 are essential for IRA2-NBD1 binding, optimal nucleotide turnover and polypeptide translocation. We propose that translocation ligands alleviate C-terminal domain suppression, allowing IRA2 to stimulate nucleotide turnover at NBD1. DEAD motors may employ similar mechanisms to translocate different enzymes along chemically unrelated biopolymers.

The Rhodobacter capsulatus genome.

The genome of Rhodobacter capsulatus has been completely sequenced. It consists of a single chromosome containing 3.5 Mb and a circular plasmid of 134 kb. This effort, started in 1992, began with a fine-structure restriction map of an overlapping set of cosmids that covered the genome. Cosmid sequencing led to a gapped genome that was filled by primer walking on the chromosome and by using lambda clones. Methods had to be developed to handle strong stops in the high GC (68%) inserts. Annotation was done with the ERGO system at Integrated Genomics, as was the reconstruction of the cell's metabolism. It was possible to recognize 3709 orfs of which functional assignments could be made with high confidence to 2392 (65%). Unusual features include the presence of numerous cryptic phage genomes embedded in the chromosome.

Genomes OnLine Database (GOLD): a monitor of genome projects world-wide.

GOLD is a comprehensive resource for accessing information related to completed and ongoing genome projects world-wide. The database currently provides information on 350 genome projects, of which 48 have been completely sequenced and their analysis published. GOLD was created in 1997 and since April 2000 it has been licensed to Integrated Genomics. The database is freely available through the URL: http://igweb.integratedgenomics.com/GOLD/.

Analysis of the Thermotoga maritima genome combining a variety of sequence similarity and genome context tools.

The proliferation of genome sequence data has led to the development of a number of tools and strategies that facilitate computational analysis. These methods include the identification of motif patterns, membership of the query sequences in family databases, metabolic pathway involvement and gene proximity. We re-examined the completely sequenced genome of Thermotoga maritima by employing the combined use of the above methods. By analyzing all 1877 proteins encoded in this genome, we identified 193 cases of conflicting annotations (10%), of which 164 are new function predictions and 29 are amendments of previously proposed assignments. These results suggest that the combined use of existing computational tools can resolve inconclusive sequence similarities and significantly improve the prediction of protein function from genome sequence.

WIT: integrated system for high-throughput genome sequence analysis and metabolic reconstruction.

The WIT (What Is There) (http://wit.mcs.anl.gov/WIT2/) system has been designed to support comparative analysis of sequenced genomes and to generate metabolic reconstructions based on chromosomal sequences and metabolic modules from the EMP/MPW family of databases. This system contains data derived from about 40 completed or nearly completed genomes. Sequence homologies, various ORF-clustering algorithms, relative gene positions on the chromosome and placement of gene products in metabolic pathways (metabolic reconstruction) can be used for the assignment of gene functions and for development of overviews of genomes within WIT. The integration of a large number of phylogenetically diverse genomes in WIT facilitates the understanding of the physiology of different organisms.

Protein interaction maps for complete genomes based on gene fusion events.

A large-scale effort to measure, detect and analyse protein-protein interactions using experimental methods is under way. These include biochemistry such as co-immunoprecipitation or crosslinking, molecular biology such as the two-hybrid system or phage display, and genetics such as unlinked noncomplementing mutant detection. Using the two-hybrid system, an international effort to analyse the complete yeast genome is in progress. Evidently, all these approaches are tedious, labour intensive and inaccurate. From a computational perspective, the question is how can we predict that two proteins interact from structure or sequence alone. Here we present a method that identifies gene-fusion events in complete genomes, solely based on sequence comparison. Because there must be selective pressure for certain genes to be fused over the course of evolution, we are able to predict functional associations of proteins. We show that 215 genes or proteins in the complete genomes of Escherichia coli, Haemophilus influenzae and Methanococcus jannaschii are involved in 64 unique fusion events. The approach is general, and can be applied even to genes of unknown function.

Universal protein families and the functional content of the last universal common ancestor.

The phylogenetic distribution of Methanococcus jannaschii proteins can provide, for the first time, an estimate of the genome content of the last common ancestor of the three domains of life. Relying on annotation and comparison with reference to the species distribution of sequence similarities results in 324 proteins forming the universal family set. This set is very well characterized and relatively small and nonredundant, containing 301 biochemical functions, of which 246 are unique. This universal function set contains mostly genes coding for energy metabolism or information processing. It appears that the Last Universal Common Ancestor was an organism with metabolic networks and genetic machinery similar to those of extant unicellular organisms.

Genomes OnLine Database (GOLD 1.0): a monitor of complete and ongoing genome projects world-wide.

UNLABELLED: GOLD (Genomes On Line Database) is a World Wide Web resource for comprehensive access to information regarding complete and ongoing genome projects around the world. AVAILABILITY: GOLD is based at the University of Illinois at Urbana-Champaign and is available at http://geta.life.uiuc.edu/ approximately nikos/genomes. html. It is also mirrored at the European Bioinformatics Institute at http://www.ebi.ac.uk/research/cgg/genomes.html. CONTACT: genomes@ebi.ac.uk

Transcription in archaea.

Using the sequences of all the known transcription-associated proteins from Bacteria and Eucarya (a total of 4,147), we have identified their homologous counterparts in the four complete archaeal genomes. Through extensive sequence comparisons, we establish the presence of 280 predicted transcription factors or transcription-associated proteins in the four archaeal genomes, of which 168 have homologs only in Bacteria, 51 have homologs only in Eucarya, and the remaining 61 have homologs in both phylogenetic domains. Although bacterial and eukaryotic transcription have very few factors in common, each exclusively shares a significantly greater number with the Archaea, especially the Bacteria. This last fact contrasts with the obvious close relationship between the archaeal and eukaryotic transcription mechanisms per se, and in particular, basic transcription initiation. We interpret these results to mean that the archaeal transcription system has retained more ancestral characteristics than have the transcription mechanisms in either of the other two domains.